Topological phonon modes and their role in dynamic instability of microtubules 
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Microtubules (MTs) are self-assembled hollow protein tubes playing important functions in live 
cells. Their building block is a protein called tubulin, which self-assemble in a particulate 2 dimen- 
sional lattice. We study the vibrational modes of this lattice and find Dirac points in the phonon 
spectrum. We discuss a splitting of the Dirac points that leads to phonon bands with non-zero Chern 
numbers, signaling the existence of topological vibrational modes localized at MTs edges, which we 
indeed observe after explicit calculations. Since these modes are robust against the large changes 
occurring at the edges during the dynamic cycle of the MTs, we can build a simple mechanical 
model to illustrate how they would participate in this phenomenon. 

PACS numbers: 63.22.-m, 87.10.-e,63.20.Pw 



MTs exist in every eukaryotic cell, being part of the 
cellular cytoskeleton and playing important roles in cell 
division and intracellular transport. MTs are made of 
a- (3 tubulin protein dimers forming long protofilaments, 
which self-assemble into hollow tubes pQ. Under constant 
chemical conditions, the tubes switch randomly between 
growing and shrinking modes [2] . This dynamic instabil- 
ity (DI) is essential to the functioning of the MTs. 

In Ref. [3 j, the MT's total energy was considered as 
function of the transversal and longitudinal curvatures of 
the tubulin sheet near the edge. It was proposed that the 
energy landscape has two minima separated by an energy 
barier, one favoring a cylindrical and the other a trumpet 
shape. Ref. [3] proposed that, during growth, the MTs 
reside near the first minimum, but somehow they start 
climbing until they overcome the energy barrier when 
they start sliding towards the second minimum, forcing 
the ends to open and the tubes to depolymerize. 

Fig. [TJa) shows a device with a similar energy land- 
scape. It is a bow made of two elastic roads hold to- 
gether by two connectors and a rubber band stretched 
between them, so that it can slide freely between the two 
roads. The energy of the bow, as function of x [Fig.[TJa)], 
has two minima. We can switch between the minima by 
pushing hard on the free end, or by using a succession of 
low energy actions. Indeed, let us attach a spring and a 
weight at the free end and aim small beads at the weight 
[Fig. []Jb)]. The energy of the falling beads is stored in 
the harmonic motion of the weight, whose oscillation am- 
plitude grows, forcing the bow at some point to climb the 
energy barrier and switch its configuration [4]. 

The spring+weight can represent a vibrational mode 
localized at the MT's edge and the beads a sequence of 
synchronized GTP hydrolyses. In ordinary lattices, the 
edge modes can easily disappear if changes occur at the 
edge. The edge mode mentioned above must be resilient, 
of special nature, because the properties of the edge vary 
wildly during DI. Switching for a moment from phonons 
to electrons, we point out that an entirely new class of 



materials, called topological insulators, has been discov- 
ered recently [5 , 6 , 7 , with the remarkable property that 
electronic states appear near any edge that is cut into 
such material. This is a consequence of the unusual bulk 
properties of the materials. For this reason, the elec- 
tronic edge states cannot be destroyed by any chemical, 
mechanical, etc., treatment of the edge [HIE]- 

Using a realistic model, we demonstrate that the bulk 
phonon spectrum of the tubulin sheet can display the 
same unusual properties seen in the topological insula- 
tors. Based on a well established connection between 
the bulk and the edge properties [8 , we conclude that 
topological phonon states can appear at the MT's edges, 
which we indeed observe after explicit calculations. We 
advance the hypothesis that these topological edge modes 
play a role in the MT's DI, similar to that of the 
spring+weight in the mechanical model of Fig. [TJa-b). 

We model the tubulin sheet as a 2D lattice of rigid 
dimers [Fig. [T][c)] with harmonic interactions [Fig.JjJd)]. 
The system can support a variety of complex motions, 
but here we want to single out a particular motion that 
displays the interesting features mentioned above. There- 
fore, we restrict this study to propagating oscillatory 
motions that involve displacements and rotations of the 
dimers in the arbitrarily chosen, but fixed xy local planes 
[Fig. flic)], with the two degrees of freedom explained in 
Fig. pTe). We assume that the degrees of freedom left 
out are involved in oscillatory motions that either occur 
in different ranges of frequencies or couple weakly with 
the oscillatory motions that we just singled out. 

The elastic energy of a spring whose ends are displaced 
as in Fig. [jjf) is ^K[e(ri-r2)] 2 [e=e/|e|], if we retain 
only the quadratic terms. Therefore, the elastic energy 
for small oscillations stored in the network of springs is: 

V = 1 /2En{Ki[e 1 (r> R+bi -r R )? 
+K 2 [e 2 (r' It+bl+h2 - r R )] 2 + K 3 [e 3 (r R+b2 - r R )] 2 
+if 4 [e 4 (r R+b2 - r R )] 2 + K 5 [e 5 (r R+h2 - r R )] 2 
+K 6 [e 6 (v' R+h2 - r R )] 2 + if 7 [e 7 (r R+b2 _ bl - r R )] 2 } . 
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TABLE I: ei,2(k)— eo are linear combinations of the first row, 
with coefficients specified in the next rows (Sj=l, except for 
j=3 and 6 when — and eo — 5^ ■ 



The displacements r (r') of the red (blue) beads can 
be written in terms of ^ and cp shown in Fig. [ljc): 
r=(^' 1 +^ 2 )i and r / =(^' 1 — £ 2 )i, where £ 2 =<i(/?. Moreover, 
since the dimers' moment of inertia is 7=2<iM, the ki- 
netic energy is: T = Er^{(4) 2 + (&) 2 }- 

The equations of motion for the Lagrangean L=T— V, 
together with the ansatz: ^ 2 = ^Re J e* k ' n-2u;t aj c ' 2 dk, 
where n= (711,712) denotes a dimer's place in the lattice, 
lead to the following equation for the normal modes: 



ei(k) -Moo 2 w(k) 
w(k) e 2 (k)-Mcj 2 



U =0, 



(i) 



where k is in the Brillouin torus [0, 2ty] x [0, 2tt]. ei,2(k) 



are given in Table [jjand [K=K(e-i 



.21. 



w(k) = (K 3 - K 6 )[l-cosk 2 ] + i(K 4 -K 5 )sink 2 
+iK\ sin k\ + zi\2 sin(fci + fe) + ^7 sin(fci — k 2 )- 



(2) 



Eq. [T] has solutions only if the determinant of the matrix 
is null and this happens only if uo takes the values: 



4(k) 



VM[e(k) ± V|Ae(k)P + Kk)|2 



(3) 



Here, e(k) and Ae(k)=|[ei(k)±€2(k)], respectively. 

Degeneracies cj_(k)=cc; + (k) occur whenever Ae(k) and 
w{k) are simultaneously zero. Due to the symmetries of 
the lattice, it is likely that Ks=Kq, in which case w(k) 
becomes purely imaginary. With this choice, we can view 
Ae(k)=0 and w(k)=0 as the equations of two closed loops 
on the Brillouin torus. There is a high probability for the 
two loops to intersect and, when it happens, degenera- 
cies appears and the two phonon bands touch each other 
at isolated points (see Fig. |2|. These special points are 
called Dirac points and their presence is a strong indica- 
tion that our model has topological properties. 

We now start the discussion of the bulk topological 
properties. If we allow a small difference between K3 
and Kq or include additional terms in the Lagrangean, 
the Dirac points split and the bands separate. In this 
case, one can define the Berry curvature: 

F±(k) = (2^)- 1 Tr{p ± (k)[S fel p ± (k),^ 2 p±(k)]}, (4) 



where p± n (k) = 



are the projectors onto the normal 



modes and [, ] denotes the usual commutator. The inte- 
grals J F±(k)d 2 k over the Brillouin torus are always inte- 
ger numbers, called Chern numbers, which remain invari- 
ant under continuous deformations of the Lagrangean, 
unless the phonon bands touch. 
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FIG. 1: The system can be switched between two equilibrium 
configurations by: (a) delivering a large amount of energy at 
once or (b) by attaching a spring+weight and releasing the 
beads one by one.(c) The 2D lattice of dimers of primitive 
vectors bi and b2. At equilibrium, the dimers are all parallel 
and orientated along an arbitrary direction. The xy coordi- 
nate system indicates the plane of motion. The y axis is along 
the dimer but the x axis is arbitrary, (d) The interaction be- 
tween the tubulin dimers is modeled by a network of springs. 
There are 7 distinct springs and corresponding e unit vectors. 
We also indicate second order neighbors whose harmonic in- 
teraction is also considered, (e) The degrees of freedom for 
the xy planar motion, (f) Stretched configuration of a spring 
with the ends displaced by 1*1,2- 



When the Chern numbers are non-zero, the system 
acquires special topological properties [10]. We show 
that this is possible in our model. As long as the La- 
grangean is symmetric under the time reversal opera- 
tion, the Chern numbers are zero. However, general 
quadratic Lagrangeans should also contain terms that in- 
volve products between velocity and displacement. Such 
terms break the symmetry mentioned above, thus allow- 
ing for non-zero Chern numbers. They can appear be- 
cause of some weak magnetic properties of the tubulin 
proteins or because of the medium surrounding the lat- 
tice. A fairly general form of such terms will be: 



y' = \Y. E E 

R m=l,2 j=l,2 

Note that a plus combination inside the accolades leads 
to a total time derivative, which can be ignored. We ex- 
pect the strength of V to be small relative to the terms 
already included in the Lagrangean. The second neigh- 
bors interactions indicated in Fig. [TJd), 

V" = V'ErE* K*[e s • (r#_ bi - r#)] 2 , (6) 
where # means primed or un-primed, could have similar 





FIG. 2: (a) The k values where Ae(k)=0 (red) and w(k)=0 
(blue), (b) Plot of the phonon frequencies u±(k). K1-7 were 
fixed at: 1.1, 0.7, 0.27, 0.75, 0.325, 0.27, 0.75, and M=l. 



strength and it should be considered at this point. 

After including V and V" into the Lagrangean, the 
equation of the normal modes takes the same form Eq. [T] 
but we must add [m=l,2] 



(£ 8 + cos + 



-\ 2 Tm,j sin kj 



(7) 



to ei 5 2, respectively, and (K 8 — i?g)[l— cos &i] to w. Here, 
c^o is a characteristic frequency that will be specified 
shortly. We now search for those values of the param- 
eters that give non-zero Chern numbers. We consider 
the parameters used to generate Fig. |2j except that now 
we allow a small mismatch Ks=0.27 and Kq=0.18. In 
addition, we take Kg =0.018 and Kg =0.1 and 71,1=- 
72,1=— 71,2=72,2=7- In contradistinction to the elec- 
tronic topological insulators, the Dirac points cannot be 
split only by the time reversal breaking term V 7 , because 
w (e(k)) remains purely imaginary (real). One should 
also be aware that the equation for the normal frequen- 
cies is also more complicated since now we have terms 
proportional to uo 2 and uo. Therefore, the observation of 
a parameter range where the Chern numbers are non-zero 
represents one of the important findings of our work. 

Let us set 7 to zero first. The Lagrangean is symmetric 
under time reversal, consequently, the Chern numbers are 
zero. We now increase 7, gradually introducing the term 
that breaks the time reversal symmetry. Even if we do so, 
the Chern numbers remain unchanged, unless the phonon 
bands touch and then separate again. We show that this 
is exactly what happens. Indeed, let us, for a moment, fix 
uj in Eq.[7]to c^o, in which case the frequencies u; T (k) are 
given by same Eq. [3j The bands touch if simultaneously: 
Ae(k)=0, Re[w(k)]=0 and lm[w(k)]=0. In Fig. ^a) we 
plot the solutions for these equations, which are closed 
loops on the Brillouin torus. We color code them with 
red, green and blue, respectively. When 7 is increased 
from to 0.5, the blue and green loops remain unchanged 
but the red loop moves in the directions indicated by the 
arrows, brushing over the intersection point ko of the 
green and blue loops. Therefore there is a 7 (=70) for 
which (j_(ko)=o;+(ko). We now unfreeze u in Eq.[7]and 
take cJo=^±(ko). This does not change the normal modes 
equation at ko, therefore the bands will still touch when 
7=70. In Figs.[3jb-c) we plot the actual phonon bands for 
three values of 7, showing how the bands touch and then 




FIG. 3: (a) The k values where Ae(k)=0 (red), Re[w(k)]=0 
(blue) and lm[iu(k)]=0 (green). The full/dash lines are for 
7=0/0.5. (b - d) The phonon spectrum for 7=0, 0.45 and 0.6. 
The last row shows the phonon spectrum of a ribbon as func- 
tion of k\\ for (e) periodic and (f) open boundary conditions. 



separate as 7 is increased. We have also verified directly 
that the Chern numbers switch their values from to ±1 
when 7 is about 0.45. 

Whenever the Chern numbers are not zero, vibrational 
modes appear near any edge cut into the 2D lattice [8]. 
We can now perform a routine phonon band calculation 
to observe these modes [7=0.6 in these calculations]. We 
considered a tubulin ribbon along the direction of b 2 and 
computed the phonon spectrum as function of ku using 
periodic and open boundary conditions (see Fig. pte-f)). 
Only in the last case one can see two distinct phonon 
bands separating from the bulk spectrum. The two bands 
are localized on the two edges of the ribbon. Also, each 
band connects the lower and upper parts of the bulk spec- 
trum, therefore, no matter how one distorts the ribbon, 
the bands will always cross the direct band gap. 

At this point we established the existence of robust 
phonon edge modes and now we return to the mecha- 
nism of DI. According to the cap theory [1 , which we 
adopt here, the bulk of the microtubules is made of pre- 
dominantly GDP-tubulin, while the edge displays a cap 
of a few rows of GTP-tubulin. The GTP and GDP 
bound tubulin prefer parallel and tilted orientations rel- 
ative to the tube's axis, respectively. During polymer- 
ization, GTP-tubulins from solution attache to the edges 
and, at almost the same rate, GTP-tubulins from the 
back rows of the cap turn into GDP-tubulins. The newly 
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formed GDP-tubulins would prefer a tilted orientation, 
but they are locked inside the lattice and therefore the 
tilt cannot occur and instead mechanical strain develops. 
One established view is that this strain energy is hold se- 
curely by the cap, but if the cap opens the strain energy 
is released and the MT depolymerize. 

All DI models are sketchy about how the MT opens 
its cap. The cap could be forced open or simply shed 
off. In either case, the MT's edge must cross an energy 
barrier as described in Ref. [3] . The only energy available 
during DI are quanta of 12 kT released whenever GTP 
hydrolysis. The energy barrier discussed in Ref. [3] was 
never measured, but a rough yet fair estimate will place 
it around 50 kT. It is then evident that the MT must 
have a way of accumulating part of the 12 kT energy and 
concentrate it near the edge. Note that the strain energy 
is uniformly distributed within the bulk of the tube. 

In our view, the sudden strain due to GTP hydrolysis 
gives a kick [TT] to the edge mode, which can enhance or 
reduce the amplitude of the mode. There is also dissi- 
pation from the edge mode into the bulk and surround- 
ing medium, therefore, we cannot expect a slow, steady 
energy buildup at the edge. Instead, the MT needs a 
small number of lucky strikes, synchronized well enough 
to push the edge mode over the energy barrier. Assum- 
ing a 50 kT barrier and a 12 kT energy transfer to the 
edge mode during the lucky strikes, the MT needs about 
5 consecutive lucky strikes to open the cap. The time 
till such strikes occur is random, which can explain the 
stochastic behavior of the catastrophe events. 

The existing calculations put the MT's vibrational 
modes in the frequency range from 1 MHz to few GHz 
and the ratio between their relaxation times and the pe- 
riod was estimate to be as high as 10 4 [T3] . The aver- 
age time between consecutive tubulin dockings can be 
derived from the polymerization speed [2] and is about 
r av =0.02 s. The kicks, however, come from random GTP 
hydrolysis at the back of the cap and involve tubulins 
already attached to the edge [12]. The average time be- 
tween consecutive hydrolyses is comparable to r av , but 
the processes themselves happen much faster and they 
are also much simpler than the docking of the tubulins. 
Therefore, the well synchronized kicks can occur within 
time intervals much smaller than r av and comparable to 
the phonon relaxation times. 

Anticancer drugs can modify DI. Taxol, for example, 
increases the growth period without strengthening the 
cap or stiffening the bulk of the MT [14]. The effect is 
present even when there is one bound taxol per several 
hundred dimers. Equally puzzling, at these small con- 
centrations, taxol stabilizes the MTs without changing 
the growth rate [15]. This means the energy flow from 



GTP hydrolysis during the growth period remains un- 
changed. It seems that the MTs edges are less efficient 
at harvesting energy in the presence of taxol. 

In our view, taxol modifies the bulk properties of the 
MTs, leading to a derealization of the edge mode. The 
localization near the edge of the topological mode is 
strong when the bulk phonon bands are widely sepa- 
rated and becomes weaker as the bands come towards 
each other [8]. In the mechanical model of Fig. [I] de- 
localization is like attaching the weight away from the 
free end of the bow. But the further we attach them 
the harder it is to switch the bow's configuration and a 
larger number of lucky strikes is needed. Thus, in our 
picture, the duration of the MT's growth depends on the 
localization of the edge mode, which is controlled by the 
MT's bulk properties. 

To conclude, we advanced the hypothesis that phonon 
edge modes play important role in the dynamical insta- 
bility of microtubules. Using an explicit lattice model, we 
demonstrated the existence of topologically robust edge 
phonon modes and we advanced the hypothesis that the 
MT uses them to concentrate energy near the edge in 
order to open the cap. The bulk lattice model can be 
extended to include an anharmonic lattice model of the 
cap, with a double- well energy landscape, allowing us to 
simulate and observe the opening of the cap explicitly. 
This and designing experiments to test our hypotheses 
are underway. 
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